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General formula describing both the divertor strike point splitting and width of magnetic islands 
created by resonant magnetic perturbations (RMPs) in a poloidally diverted tokamak equilibrium 
is derived. Under the assumption that the RMP is produced by coils at the low-field side such eis 
those used to control edge localized modes (ELMs) it is demonstrated that the width of islands 
on different magnetic surfaces at the edge and the amount of divertor splitting are related to each 
other. Explanation is provided of aligned maxima of the perturbation spectra with the safety factor 
profile - an effect empirically observed in models of many perturbation coil designs. 
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I. INTRODUCTION 

Resonant magnetic perturbations (RMPs) are being 
investigated as tool to control edge localized modes 
(ELMs), in particular their application to ITER is fore- 
seen. The RMPs for ELM control are produced by 
coils whose design is specific to each tokamak. This 
method is reminiscent of the ergodic divertor on toka- 
maks Tore Supra and TEXTOR which also essentialy 
relies on RMPs produced by dedicated systems of coils, 
and many aspects are similar, especially the formation 
of magnetic islands on rational surfaces and possible 
stochastic transport when the islands overlap. However 
the presence of X-point in the poloidally diverted toka- 
maks provides some unique features: the splitting of di- 
vertor strike points, or divertor footprints ~ a signature 
of the homoclinic tangle created by the perturbation-'' , 
and the divergence of the safety factor profile at the sep- 
aratrix, due to which the number of rational surfaces 
even for one toroidal mode is infinite and island over- 
lap on them is facilitated. The noncircular geometry 
of the plasma cross-section also complicates analytical 
treatment of magnetic islands, requiring cautious use of 
non-orthogonal coordinate systems. 

Since both the homoclinic tangle and magnetic islands 
with the resulting stochasticity are consequences of the 
perturbation, it is natural to ask if there is any relation 
between them. It has been known that while every island 
chain is related to a different polodal mode of the pertur- 
bation, their sizes can be expressed by a single function 
- the Poincare-type integral which also generalizes to an 
expression for the length of the divertor footprinta^.The 
present paper uses instead the more familiar formalism of 
the Melnikov integral^ (a standard tool for the analysis 
of the homoclinic tangle) and explores the radial depen- 
dence of this function in the often encountered case of 
RMPs from the perturbation coils localized at the low- 
field side (LPS), as it is the case for most ELM control 
coil designs. 

The paper starts with a review of the method to calcu- 
late island sizes. Care is taken to provide a formula valid 
in a plasma with a general aspect ratio and non-circular 



cross-section as this is essential for the edge pedestal re- 
gion near the separatrix which is crucial for ELM con- 
trol. It is shown how the island sizes are determined by 
a function whose definition is not affected by the coor- 
dinate singularity on the separatrix, and which we call 
the Melnikov-like function. The homoclinic tangle and 
divertor footprints are then explained together with the 
method of Melnikov integral for determining analytically 
the divertor footprints length. A particularily simple ex- 
pression is given for the case when the perurbation has 
only one dominant toroidal mode. The relation between 
Melnikov integral and the Melnikov-like function is ex- 
plained. Then we restrict our treatment to the case of 
perturbation localized on the LFS. Under this assump- 
tion the relation between modes at the different surfaces 
is derived and it is shown how does the divertor footprint 
length relate to the sizes of magnetic islands at the edge. 



II. WIDTH OF MAGNETIC ISLANDS IN A 
REALISTIC GEOMETRY 

A. Expression using the resonant modes of the 
perturbation 

Nonaxisymmetric magnetic perturbations of a toka- 
mak magnetic field produce chains of magnetic islands 
on magnetic surfaces with low-order rational values of 
the safety factor q. Those islands are created by reso- 
nances of the perturbation field ^B with the unperturbed 
field lines on the rational surface. To express these res- 
onances we use a magnetic coordinate systenP (s, 0*,^) 
where s is a flux surface label, 9* the poloidal coordinate 
and if the toroidal coordinate. The angular coordinates 
9* , ip are chosen as in the PEST coordinates^ and the ra- 
dial coordinate s is defined as square root of the normal- 
ized poloidal magnetic flux tpN: s = \/0JVj in accordance 
with previous works (e.g. Refs. I6l7p . The coordinate (p 
is simply the geometric toroidal angle and the coordi- 
nate system has the property that 9* = oip/q{s) + const, 
along a fleld line, where o represents the orientation of 
the magnetic fleld: o = 1 for the left-handed fleld and 



o = — 1 for the right-handed fieltr^. The islands are cre- 
ated by the (m, n) Fourier component b}^ ^^ of the nor- 

mahzed perturbation Sb^ = B^ /B^ , where B^ = (5B • Vs 
is the contravariant s component of the perturbation and 
B"^ = B • V(/9 is the contravariant (^ component of the 
equihbrium field. The Fourier transform is taken with 
respect to the 9* and ip coordinates, thus we have 



5h^ 



E 



^(m,n) 



exp [i{m9* — riLp)] 



(1) 
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and the Fourier harmonics can be obtained as 



hUn) = (^ / exp[-i(mr - n^)W d^ d^. (2) 

A Fourier component is resonant with the unperturbed 
field lines when q — om/n. The values b}_^ ^ n are com- 
plex, and as 5b^ is real, the following relation for the 
complex conjugate holds: (&(„^„)j = &(_„_„)■ An 

alternative to ([I]) is a representation using purely real 
coefficients: 



(jn.n) 



Sb^ = Y^ 6[„_„) sin(TO6'* - nip + Xmn) (3) 

m— — oo,n— 1 
2|^(m,n)l 



(4) 
(5) 



Widely used formulae exist for determining the width 
of magnetic islands from the Fourier spectrum of the 
perturbatiorPEl. They are typically derived in a cylin- 
drical geometry where the toroidal curvature is not be- 
ing taken into account (the toroidal magnetic field is 
considered constant), thus B^ in the expression for Sb^ 
is approximated by its value at the magnetic axis. As 
noted in Ref "9", this leads to an error in estimating the 
island width. In the example of TEXTOR and its dy- 
namic ergodic divertor (DED), the island size was over- 
estimated because the DED coils are located at the high- 
field side, thus the actual value of B^ is larger and Sb^ 
is smaller than in the cylindrical approximation. For the 
ergodic divertor of Tore Supra, which was located at the 
low-field side, the island sizes were underestimated. It 
should be noted that while the toroidal field magnitude 
_Bt varies with the radial distance R from the major axis 
as -Bt oc 1/-R, for the contravariant component the de- 
pendence is stronger: B^ ex 1/B?. 

Moreover the cylindrical formula for island widths uses 
r (the distance from the magnetic axis) as a radial coordi- 
nate and thus is valid only in a situation where the mag- 
netic surfaces have circular and concentric cross-sections. 
In divertor tokamaks we are far from this geometry, espe- 
cially in the edge region near the separatrix which is the 
most important when perturbations are used as an ELM 
control mechanism. We thus need a formula which would 
be usable in a general geometry, with a varying toroidal 
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field and noncircular flux surfaces, using for example the 
coordinate s as a general flux surface label instead of r. 
To derive this formula we introduce new coordinates 
X = 0* — n/mp and s = s — sq where sq is the flux label 
of the resonant surface where q = m/n. The differential 
equation of the field line are 

Dl 

(6) 

(7) 

Using the coordinates x ^^^d s and the relation 

dx ^ 1 _ n (3^ 

dp q m 

the equation ([6]) becomes 

l/q - n/m ds ^ B^/B^ dx (9) 

Keeping only the resonant part of the perturbation, thus 
substituting B^/B^ by 6j^^^ sin(mx + Xmn), we obtain 

l/q - n/m ds = 6(„ „) sin(mx + Xmn) dx (10) 

Using a linear approximation of the left side, we obtain 
dq^ 



ds 



s ds w 6(„^„) sin(mx + Xmn) dx (H) 

q—m/n 

This equation can be easily integrated to obtain an alge- 
braic equation for field lines: 

. 2q'b^ 



[cOs(toX + Xmn) + C] 



(12) 



(rn.n) 

s- w ^^ 

q'm 

where q' = dq/ ds at the resonant surface and C is an 
integration constant. The choice C = 1 corresponds to 
the island separatrix whose maximum radial excursion is 
the island half-width i5, given by the formula 



S = 



/Vfei 



(m,n) 



q'm 



(13) 



An alternative to this approach is to use a Hamilto- 
nian approach where the field lines are interpreted as 
trajectories of a Hamiltonian dynamical system whose 
Hamilton function is the poloidal flux and the pertur- 
bation is represented as a perturbed Hamiltonian (flux). 
This approach has been used in many theoretical works. 
We briefly review it in the appendix [X] and prove its 
equivalence to the approach described above. It should 
be emphasized that the hamiltonian approach automat- 
ically includes correctly the effects of toroidal geometry 
and non-circular cross-section - no corrections are nec- 
essary. It is however still important to have a correct 



formula using the perturbed magnetic fleld ( 13 ), because 



this is the approach usually used in numerical studies of 
perturbation coil designs, as the perturbed field can be 
readily calculated from the coil geometry by the Biot- 
Savart formula. We will also see that the harmonics of 
the perturbed field are directly related to the Melnikov 
function. 



B. Expression using Melnikov-type integral 

The coordinate system (6**, s) on a poloidal plane has a 
singularity at the separatrix. It is useful to define a value 
characterizing magnetic islands which, unlike bl s , will 
not use the 9* coordinate, so it will stay well-defined even 
at the separatrix. 

We start by defining a coordinate (j) which will be used 
instead of 9* . We follow a procedure used in the defini- 
tion of the separatrix map and the Melnikov integraP'^^. 
For every magnetic surface the point on the outboard 
midplane has (f> — 0. Following a field line parameter- 
ized by the toroidal angle ip from this point, we assign 
to any other point on the field line the value (j) = ip. 
Thus (p oi a given point is the toroidal angle needed to 
reach it by following a field line from the outboard mid- 
plane. Since the field line returns to the same poloidal 
position after making q toroidal turns, the range of the 
coordinate (p needed to cover a magnetic surface in the 
poloidal plane is (—qn^qn) where the endpoints of this 
interval are identified with each other. Together with a 
flux surface label such as s we obtain a coordinate system 
on the poloidal plane. The separatrix is a special case: 
it is covered hy (j) € (— oo, c») since q is infinite on the 
separatrix, and the X-point corresponds to ^ = ±cx). In 
this case (p is called a homoclinic coordinate^. As on a 
field line ip = oq9* + const, and 6** = on the outboard 
midplane where (f) = 0, the relation between 9* and (j) is 
(f> = oq9* . Using this relation, the definition i2n of 6J . 
can be rewritten as 



bl 



(m.n) 

1 

oq{2ny 



q-K /'27r 
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exp — 1 [ o — (p ~ nip I 5b dipAcj) 

(14) 
If we define a toroidal perturbation Fourier mode b\ as 

~b'ni^)^^J eMinip)6b\^,cp)dv (15) 



we may write ( 14 ) as 
1 /"'"" 



bli^)d<P (16) 



On resonant surfaces with q = om/n this may be simpli- 
fied to 



(m^onq,n) 



1 



o2TTq J_y^ 



exp(-in0)6,\(0)d0 (17) 



{b\, bl . and 6b^ all depend also on the magnetic sur- 
face, this was omitted from the expressions above for 



brevity). In (17) there appears a complex Melnikov-like 
function S'„(s) given by 



qTT 



Sn{s) = I exp(-i?i0)6^(s,( 

-qTT 



(18) 



defined using the coordinate (p which does not have a 
singularity at the separatrix, so the definition can be ex- 
tended to the separatrix: 



/oo 
cxp{-mcj))bl,{s = l,^)dcf,. 
OO 



(19) 



The function 5'„ fulfills our requirement: it can replace 
^}m n) ~ Sn/{o2'Kq) and is defined using values which 
remain regular at the separatrix. The island width ( 13 ) 
can be expressed using Sn instead of bl ^^ : 



6 = 



i Ms„\ 

mrq' 



(20) 



The only remaining divergent term in ( 20 ) is the shear 



q' which grows to infinity at the separatrix. This depen- 
dence is physical: its consequence is that island width 
has a zero limit at the separatrix. 



III. DIVERTOR FOOTPRINTS 

Since the particle and heat transport are mostly paral- 
lel to field lines, the patterns of particle and heat flux to 
the divertor plates can be expected to be related to the 
divertor magnetic footprints, i.e. the patterns of inter- 
sections of field lines with the divertor. Field lines which 
carry heat and particle fluxes from inside the plasma are 
those with a high connection length, i.e. the number of 
toroidal turns following the field line in the plasma before 
it reaches the wall again. 

Since the field lines can be interpreted as trajectories of 
a Hamiltonian dynamical system with the toroidal angle 
in the role of the time, methods of the theory of Hamil- 
tonian systems can be used. A concept especially useful 
for the study of divertor footprints is the one of invari- 
ant manifolda-3. An invariant manifold is a surface in the 
phase space of the dynamical system which remains in- 
variant by the time evolution of the system, thus a trajec- 
tory with an initial point on the invariant manifold is con- 
strained to remain on iiPl. In our case the trajectories are 
field lines and one example of invariant manifolds are the 
magnetic surfaces of the toroidally symmetric tokamak 
equilibrium. A particularly interesting case of invariant 
manifolds are the stable and unstable manifolds of hyper- 
bolic fixed points. A stable manifold is formed by field 
lines asymptotically approaching the fixed point, while 
the unstable manifold is formed by field lines asymptot- 
ically leaving the fixed point. The definition depends on 
the direction in which the field lines are followed. If we 
follow them in the opposite direction, the stable mani- 
fold becomes unstable and vice versa. In plasma equi- 
libria the hyperbolic fixed points are called X-points and 
are associated with the poloidal divertor or with mag- 
netic islands. An example of invariant manifolds to a 
fixed point is the separatrix of a toroidally symmetric 
configuration with a poloidal divertor. Here the stable 



and unstable manifolds coincide to form the separatrix. 
When a perturbation appears, the separatrix splits into 
the stable and unstable manifolds which no longer co- 
incide, but intersect transversally infinitely many times. 
Close to the X-point in the direction from which the field 
lines approach it (the stable direction) the stable mani- 
fold is close to the unperturbed separatrix, but the unsta- 
ble manifold wildly oscillates, creating lobes that become 
longer and narrower when the X-point is approachecP. In 
the direction of field lines leaving the X-point (the unsta- 
ble direction) we obtain a similar picture with the roles 
of the stable and unstable manifolds reversed. This com- 
plex structure is called a homoclinic or heteroclinic tan- 
gle. An important property of the invariant manifolds is 
that field lines can't cross them, because field lines can't 
intersect. Invariant manifolds thus act as boundaries for 
the field lines. Field lines originating in the hot plasma 
core are contained inside the invariant manifolds of the 
X-point and the only way they can reach the divertor tar- 
gets is when the lobes of the (un)stable manifolds near 
the X-point intersect the target plates^. By tracing the 
intersection of the manifolds with the plates one obtains 
curves which delimit the region connected to the plasma 
core, characterized by mostly high connection length in 
the laminar plot. Those divertor footprints typically take 
the form of long spiralling bands, each band correspond- 
ing to the intersection of a protruding lobe of a stable or 
unstable manifold with the divertor. 

The divertor footprints have a complicated inner struc- 
ture and not all points inside the manifolds have high 
connection length. Some of them are connected to the 
opposite divertor plate after two poloidal turns by lami- 
nar flux tubes which do not penetrate deeply under the 
separatrix. The points with high connection lengths are 
the images of invariant manifolds of the X-points of the 
inner island chains^ ^. This fine structure was studied in 
detail irP^. Here we focus on the on the overall shape 
of the divertor footprints which is given by the invariant 
manifolds of the divertor X-point and is better experi- 
mentally accessible. 

The length of the spiral can be characterized by the 
maximum value of s reached, i.e. the value stip at its 
tip. The difference Asmax of Stip and the separatrix 
value s = 1 expresses the radial distance on the divertor 
plate between the footprint's tip and the unperturbed 
strike point, which lies at the intersection of the unper- 
turbed separatrix with the divertor plate. The unstable 
manifold^ is the footprint's boundary and so the foot- 
print's tip is the point on the manifold which is the most 
distant from the unperturbed separatrix, the distance in 
terms of s being Asmax- The value Asmax thus quantifies 
the magnitude of the separatrix splitting. 

To estimate the separation between the unperturbed 
separatrix and the unstable manifold we will follow two 
field lines - one in the unperturbed field, lying on the 
separatrix, and the other in the perturbed field, lying 
on the unstable manifold. They are parameterized by 
the toroidal angle (p. Let us choose them so that they 



are initially (in the vicinity of the X-point which they 
approach asymptotically when followed backwards in ip) 
close to each other. The parametric equations of the 
perturbed field line are 






(21) 
(22) 



For the unperturbed field line they can be written ex- 
plicitely using the definition of (j) and the fact that the 
unperturbed field line lies on the separatrix where s = 1: 



s = s{ip) — 1 



(23) 
(24) 



where the constant (f>o determines the toroidal phase: it is 
the toroidal angle of the point where the field line crosses 
the outboard midplane. 

The rate of change of s along the perturbed field line 
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(25) 



If the perturbed field line does not deviate significantly 
as a result of the perturbation, we may use a first-order 
perturbative approximation and evaluate b^ on the un- 
perturbed field line: 



ds^ 

dp 



6bHs = l,cj){p,),p) 



(26) 



The deviation As of the perturbed field line from the 
unperturbed one after a full poloidal turn is given by the 
integral of (26): 



As(</.o)=/ Sb\s = l,^p)=ip-cl)o.p)dip (27) 
or using (p as the parameter: 

/•ex: 

As{cf>o) = SiM = 6b\s = l,<l>,^ + cf>o)dcb. (28) 



Note that As((/)o) is a function of the toroidal phase (po 
and may be zero: this happens when the unstable mani- 
fold intersects the unperturbed separatrix. 



The function S defined by the integral (28) is closely 
related to the Melnikov function M: 



MiM 



Sh^'ii 



dip 
d^ 



S{M (29) 



where 56'^(s = 1, </), + (/)o) = ^h\s = l,0,0 + 0o) is 
the contravariant component of 5B with respect to the 
-0 coordinate. The only difference between M and S is 
that S gives the change of s while M gives the change of 



If there is only one toroidal mode &^((/>) of the pertur- 
bation [cf. equation (15l], the function S can be replaced 



by a single complex number Sn{s = 1): 

/•OO 
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2^ |exp[-in(0 + (t)o)]bl{s = 1,0)} d0 
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exp(— in(/)o) / exp(— in(/))6,\(s = 1, 



exp(-in0o)<S'n(s = 1) 



(30) 



Sn{s = 1) is defined by the equation (19), which also 
naturally extends the definition of S'„(s) to the domain 
s < 1. Analogously the value Af„(s) = ^<Sn(s) can 
be used to express the Melnikov function M as a single 
number Mn{s = 1). 

The value Asmax is the maximum deviation of the un- 
stable manifold: Asmax = rnax^^ As{(I>q). Using S'„(s = 
1) it is expressed as 



2|^„(.s = l)| 



(31) 



The island widths and the magnitude of separatrix split- 
ting, and consequently the length of the divertor foot- 
prints, are given by a single function S'„(s): the island 
widths by its values at s < 1 and the magnitude of split- 
ting by the value at s = 1. 

It has been already known that island widths and the 
magnitude of splitting can be described by a signle ra- 
dial function^ : the Poincare-type integral i?„ which is an 
integ ral o f the modes of the perturbed poloidal flux Hi 
(cf. (A8)) instead of the perturbed field 6^. It can be 



shown that Mn — enRn so for a single toroidal mode of 
the perturbation our formalism of the Melnikov- like func- 
tions 5* or M is equivalent to the Poincare-type integral 
approach. The reference [2] gives also other results ex- 
pressed in terms of the function i?„ such as the width of 
the stochastic layer and the field line diffusion coefficient 
which can be also simply reformulated using Melnikov- 
like functions. 



IV. SPECIFIC FORM OF THE MODES OF A 
LOCALIZED LOW-FIELD SIDE PERTURBATION 

In the previous sections we introduced the Melnikov- 
like function S'„ and showed how it expresses both the 
sizes of the magnetic islands and the sizes of the diver- 
tor footprints. This is not sufficient to relate the sizes of 
the footprints to the sizes of the islands unless the radial 
dependence of Sn is known. In this section an approxi- 
mative form of this dependence at the edge will be given 
for the special case of external magnetic perturbations 
imposed by coils located at the low-field side. The mo- 
tivation for this case is the use of such coils as an ELM 



control mechanisrrP^ where the the coils are supposed to 
impact the edge region where the ELMs originate. 

We will use a simplified model of the perturbed mag- 
netic field where the perturbation is localized at the low 
field side where the field line pitch angle dip/ d6 {6 being 
the geometric poloidal angle) is assumed to be constant 
poloidally and radially. This is a realistic assumption 
for the edge region near the separatrix which is our re- 
gion of interest. We will note the local pitch angle gi: 
qi = dip/ dO = const. The variation of the safety factor is 
assumed to be caused only by the variation of the pitch 
angle in the regions where the perturbation is negligible: 
the high-field side and especially the X-point. This re- 
quires the perturbation coils to be placed sufficiently far 
from the X-point region. 

Along a field line in the low field side region we have 

e* = oip/q ^ oqiO/q. (32) 

It follows that the <j) function has a simple dependence 
on 9 in this region: 



(jj^qii 



(33) 



The m Fourier component of the perturbation w.r.t the 
geometric poloidal angle 9 is defined as 






1 

2^ 



exp{~ifn9)b'^{9) dO 



(34) 



where b'^{9) is the n toroidal Fourier component of 
6b^{9, if) considered as a function of 9: 



b'^{9)^bl{<l>{9)). 



(35) 



We will now find the relation between the Fourier com- 



ponents &[,„_„) and 6; 



(m,n) 



We are neglecting the pertur- 



bation outside the region where the Eq. ( 33 ) holds which 
allows to express bl . [Eq. ( 16 )] in terms of the 9 coor- 
dinate: 
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o2TTq J_^ 



exp 
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bni0)d9 (36) 



Equations ( 36 1 and ( 35 ) finally give simple relations be- 



tween 6i„^_„) and 6'(Lj „) 



11 = ^ft'i = '^h'^ 

(™:») q (m,qi/q,n) "(ngi.n) 



(37) 



and 



(m.n) (raqjqx.n)- 

From those it can be seen why the maxima and min- 
ima of the spectrum 6J ^n in (m, s) space form "ridges" 
and "valleys" aligned with the q profile, as can be seen 
e.g. for the proposed ITER designs in Ref. [7] (Fig. 15c 
therein) and noted for DIII-D in Ref. [M] (see Fig. lb 
therein). The Fourier component w.r.t. 9* - b}^^-. - 



is given by Eq. (37 1. Assuming that the Fourier com- 
ponent of the perturbation w.r.t. 9 - b',^ ^s - does not 
change significantly between different magnetic surfaces, 
the only radial dependence is the inverse proportional- 
ity to q which is the same for all the poloidal modes. If 
^}m n) ^^^ ^^^ maximum on one surface with s — si for 
m = mmax(si), on other resonant surface with s — S2 
it will have also maximum for m = TOniax('S2) equal to 
'>T^nia.xisi)qis2) / q{si) SO maxima will be aligned with the 
q profile which is given in the (m, s) space as the set of 
points satisfying m = nq{s). 

Using these resuls the an approximate radial depen- 
dence of bl J^^ can be found. The radial dependence 

of the geometric poloidal Fourier component B^ of the 
radial perturbation B'' = SH ■ e,. (with e^ being the 
unit vector perpendicular to the magnetic surfaces) is^S 
B^ ex r™~^. The contravariant s component B^ is given 
by B^ — B^ds/dr. Assuming that ds/dr and B^ do 
not depend significantly on the poloidal angle in the 
area with a non-negligible perturbation, the geometric 
poloidal Fourier component of the normalized contravari- 
ant perturbation is given by 



(rn,n) 



CX r 



drB^' 



(38) 



The radial dependence of b] -, is given by the for- 



{m,n) 



mula (37) where Eq. (38) can be used to substitute for 

{mqi/q,n)' 

For resonant modes we may use the Melnikov-like func- 
tion instead. From Eqs. ( 33 ) , ( 18 1 and ( 34 1 it follows that 



'^"^''^ " oi^("«i^") 



(39) 



The radial dependence of Sn{s) can be obtained from 
Eqs. (ISSl) and (pgl): 



~Sn{s) cx gir'"?!- 



IhB^ 



(40) 



At a sufficiently narrow edge region the right-hand side 
is not strongly radially dependent, so we may expect 
the values of Sn{s) on different resonant surface to be 
strongly correlated. Note that (40) and this conclusion 
applies also to the value on separatrix 5„(s = 1), which 
is thus the limit of Sn{s) at the resonant surfaces ap- 
proaching the separatrix, because ( 40 ) does not contain 
discontinuous terms. 



it especially suitable for estimating the edge ergodization 
in an X-point tokamak geometry, where the edge region is 
substantially different from a cylindrical approximation. 
The formula is based on the perturbed magnetic field and 
we demonstrated its equivalence to formulae expressed in 
terms of the perturbed poloidal fiux. We then formulated 
assumptions about the form of the perturbed magnetic 
field which correspond to the perturbations typically used 
inmo the ongoing effort to control ELMs with magnetic 
perturbations on a range of tokamaks. Namely, we sup- 
pose that the perturbation acts mostly in a region away 
from the X-point, where the pitch angle of the field lines 
does not have a significant radial variation in the region of 
interest, which is the edge zone near the separatrix. This 
assumption is valid for the coils used for ELM control 
experiments in most tokamaks, as well as the proposed 
coils for ITER. Using this assumption we then derived 
more concrete results about the perturbation harmon- 
ics which determine the island sizes. We demonstrated 
that all the resonant harmonics are correlated. Our re- 
sult expresses formally the alignment of the maxima and 
minima of the perturbation spectra with the safety fac- 
tor profile, which is often observed in the calculations of 
perturbation harmonics. We also show that the quantity 
which determines the island sizes is also directly linked 
to the Melnikov integral and thus determines the extent 
of the footprints on the divertor plates. 

Our results show that by using coils on low field side 
it is not possible to create significantly diffferent reso- 
nant perturbations on different rational surfaces. Max- 
imizing the resonant mode on one surface also leads to 
maximization of resonant modes on other surfaces. This 
is advantageous if one wants to optimize the coil sys- 
tem for maximum island overlap and stochastization. If 
one rather wants to study the effect on perturbation on 
each surface separately it might be more advantageous 
to choose a different position of the coils, as it is the case 
for the new perturbation coils on DIII-D. Maximizing the 
island overlap will also lead to maximization of divertor 
footprints due to the relation between island sizes and 
the Melnikov integral. 

As our method is restricted to a LFS-localized pertur- 
bation, the results do not apply to a perturbation field 
created inside the plasma itself, e.g. a locked mode. In 
this case the relation between magnetic islands and the 
divertor footprints may be much less constrained. 
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Appendix A: Hamiltonian representation of field 
lines and magnetic islands 

In the theory of hamiltonian dynamical systems (see 
Ref . ITB|) ■ the formula for the island width is derived using 
the hamiltonian description of field line dynamics, with 
the poloidal flux function in the role of the hamiltonian 
and the toroidal angle in the role of time (see e.g. Ref. [2]). 
The hamiltonian is defined as 



H = A., 



RA^ 



(Al) 



where A^ is the covariant toroidal component of the vec- 
tor potential and A(p = A • e^ is the physical component, 
with e^ being the unit basis vector in the toroidal di- 
rection. A convenient choice of canonical coordinates is 
the action-angle representation, where the action is the 
toroidal flux <& and the angle is 9*. The Hamiltonian 
equations are: 



AO* _ dH 

d^ ~ a¥ 

d$ _ dH 

d^ ~^de* 



(A2) 
(A3) 



In the equilibrium case H is a function of poloidal po- 
sition only and is independent on the toroidal angle (p. 
Moreover, $ and 9* being action-angle variables, they are 
chosen so that H is only a function of $ and Eq. (A3) 
is identically zero. A nonaxisymmetric perturbation is 
represented by the addition of a small term eHi ($, 0* , (f) 
to the hamiltonian, which can then be written as 



iJ($,r,(^) = Hoi'P)+eHi{'S?,d*,ip). 



(A4) 



Hq is the equilibrium part, which can be identified 
with the unperturbed poloidal flux ip. The per- 
turbed part eHi{^,0* ,Lp) corresponds to a perturbation 
dA^{^, 9* , if) of Aip. The equilibrium part has the prop- 
erty 



dHo 
d$ 



dV^ 
d¥ 



which reduces Eq. ( A2 ) to the form ajj . (We assume 
that the perturbation term e — ^^^r^ — — is negligible in 



comparison with the equilibrium term 1/q and can be 
neglected.) The equation (6]) can be derived from (A3) 
by expressing the perturbed field Sb^ using the perturbed 
potential SAip. This expression is 



„, 1 _ ds 1 86 A^ 
dip q do* 



(A6) 



The derivative ^ can be expressed as 

ds ds dV^ d$ ds 1 d$ 

dkp dib d $ di^ dip q dtp 

From Eqs. (A6| and (A7) it follows that Eq. ^ is equiv- 
alent to ( A3 1 . 



(A7) 



It is useful to decompose the perturbed potential in 
Fourier modes, analogously to the decomposition ([I]) of 
Sb^: 



dA^ 



eHi = e ^ -ff(m,„) exp [i{m6* - n</3](]A8) 



m.n— — OO 



e^ -ff(m,n) COs{m0* -nip + Xmn) 



(A9) 



From (IT]), ( |A8[ ) and (A6) we obtain the relation between 
^(m,„) andi/(„,„): 



bl 



{m,n) 



b). 



ds 1. 

d-ip q 
ds 1 



-vmeH, 



{■m^n) 



{m,n) 



dip q 



-meH, 



[m^n) ■ 



(AlO) 
(All) 



The half- width of islands measured in terms of the ac- 
tion variable (toroidal flux $) is^^. 



2q, 



leH, 



{m,n) 



dq 



(A12) 



(A5) ji^ a linear approximation, the half-width in terms of s 



is related to (5$ by the relation Sg = ^<5$. Moreover, 



dq 
d* 



d$ 



q' and 



d$ 



j4-, so using Eq. (|All|) we see 



that the expressions (A12) and (13) are equivalent. 
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